Species Numbers versus Area in Chowdhury Ecosystems 
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Abstract: For a fixed number of food levels on a lattice, the Chowdhury 
model for ecosystem is simulated and found to give a number of different 
. I species which approaches exponentially a plateau if the considered area is 

.2 ' increased. 

The Chowdhury model PP is one of several individual-based models for 
ecosystems with many (prey-predator) species |21 E] • We use it here with the 
parameters of jlj, that means a low speciation rate 0.001, and low diffusivity 
0.001, and a fixed number (usually six) of trophic layers. We refer to [TIIHIH] 
for more details on this model. Here we check if the number of different 
species which are found in a square of size R x R increases with a small 
power of the area R^, as summarized in 

So after some equilibration time (one sixth of the total number of iter- 
ations) we sum at each iteration over all different species found in a corner 
of size Rx R of the L X L lattice, with 1 < i? < L and 10 < L < 100. 
To see better the trends with increasing R we look at the i?-derivative of 
^1 this number of species, that means we determine how many new species were 

^ ' found in a square of side length R which were not already found in its inner 

^ ■ square of side length R — 1. A species in this asexual model is defined to 

consist of all animals occupying the corresponding ecological niches on the 
considered lattice sites. 

Figure 1 shows that this number of new species decays roughly exponen- 
tially with increasing R (though for small R also a power law fits). Thus the 
total number of species approaches exponentially a plateau for large R, in 
contrast to the desired power law 6J. This exp(— i?/^) decay with R is found 
to be independent of the simulation time t, Fig. 2a, while the decay length 
^ increases when the number of food levels increases from 3 to 7, Fig.2b. 
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Figure 1: Number of new species in i?-squares which are not yet found in 
the inner squares of size R — 1, summed over 20,000 time steps (iterations) 
for lattices of hnear dimension L = 10 (+), 20 (x), 50 (stars), and 100 
(squares) from bottom to top (one sample each, except sum over ten samples 
for L — 10). The counted numbers are divided by the number 2R — 1 
of involved lattice sites and rounded down to integer values. Each lattice 
carried six food levels with 1 to 32 ecological niches. 

An inhomogeneous ecosystem with different but fixed numbers of food levels 
thus might give a power law as a superposition of many exponential laws. 
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Figure 2: As Fig.l but in part a for t = 2,400 to 2,400,000 iterations at 
L = 50, and in part b for three to seven food levels (from bottom to top) at 
L = 50, t = 24,000. 
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